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Abstract 

The three-gluon glueball states are studied with the generalization of a semirelativistic potential 
model giving good results for two-gluon glueballs. The Hamiltonian depends only on 3 parameters 
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Qh' fixed on two-gluon glueball spectra: the strong coupling constant, the string tension, and a gluon 

, ' size which removes singularities in the potential. The Casimir scaling determines the structure of 
the confinement. Our results are in good agreement with other approaches and lattice calculation 
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I. INTRODUCTION 



The Quantum Chromodynamics (QCD) theory allows the existence of bound states of 
gluons, called glueballs, but no firm experimental discovery of such states has been obtained 
yet. Glueballs are important in the understanding of several mechanisms. Some authors 
proposed that glueballs characterize de-confinement [l| and that their coupling to proton 
affects the gluonic contribution to proton spin 

An important difficulty is that glueball states might possibly mix strongly with nearby 
meson states Nevertheless, the computation of pure gluon glueballs remains an inter- 
esting task. This could guide experimental searches and provide some calibration for more 
realistic models of glueballs. 

Lattice calculations are undoubtedly a powerful tool to investigate the structure of glue- 
balls. A previous study J4] predicts the existence of a lot of resonances between 2 and 
4 GeV. A recent update of this work 5| confirms the results already obtained. It is worth 
mentioning that recent lattice calculations confirm the hierarchy of the glueball spectrum 
fj. 

The potential model, which is so successful to describe bound states of quarks, is also a 
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possible approach to study glueballs [7J, |8|, |9 
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13|. In a recent paper [14j, a semirel- 



ativistic Hamiltonian is used to compute two-gluon glueballs with masses in good agreement 
with those obtained by lattice calculations of Ref. J4J1. This Hamiltonian, the model III in 
Ref. [ijj], relies on the auxiliary field formalism 15, 16] and on a one-gluon exchange (OGE) 
interaction proposed in Ref. 7J]. It depends only on three parameters: the strong coupling 
constant as, the string tension a, and a gluon size 7 which removes singularities in the 
short-range part of the potential. The constituent gluon mass is dynamically generated and 
it is assumed that the Casimir scaling determines the color structure of the confinement. 
These two ingredients are actually necessary to obtain a good agreement between the results 
from a potential model and from lattice calculations. 

In a previous paper 17], we generalized the model built for two-gluon systems in Ref. [ijj] 
for the low-lying (L = 0) spectrum three-gluon glueballs. The purpose of this paper is to 
extend the results presented in Ref. [l7| to higher glueballs. Compared to previous models 
, our approach is characterized by some improved features: semirelativistic kinematics, 
more realistic confinement, dynamical definition of the gluon mass, consistent treatment of 
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the gluon size. These points will be detailed below. The masses of the lowest glueballs 

n n n 

are computed with a great accuracy and compared with lattice calculations [4J, [5J, [6J. In 
Sec. HIl the three-gluon Hamiltonian is built, and the structure of the studied glueballs is 
presented in Sec. IIHI The three-gluon glueball spectrum is presented with the two-gluon 
glueball spectrum from Ref. Jjj] and is discussed in Sec. IIVI Some concluding remarks are 
given in Sec. EJ 



II. HAMILTONIAN 



A. Parameters 



In Ref 



two sets of parameters, denoted A and B, were presented for the model III 



(see Tabled]). With the set A, it is possible to obtain glueball masses in agreement with the 
results of some experimental works Q, 19]: the lowest 2 ++ state near 2 GeV, the lowest 



++ state near 1.5 GeV, and the lowest + state near 2.1 GeV. The values of a and as 
are close to the ones used in some recent baryon calculations 



2I|. With the set B, glueball 

masses were computed in agreement with the results of the lattice calculations of Ref. {4]. 



If the absolute glueball masses found in Ref. 14J with both sets are strongly different, the 
relative spectra are nearly identical. As we use in this work a three-body generalization of 
the Hamiltonian model III of Ref. [141 ]. the two sets will also be considered. Nevertheless, 
in the following, we mainly focus our attention to the results obtained with the set B. 

It is worth mentioning how the parameters have been determined in Ref. jl4J. The mass 
of the lightest 2 ++ is nearly independent of the values of 0:5 and 7, but depends strongly 
on a. So, this last parameter has been determined with this 2 ++ state. The remaining 
parameters as and 7 have then be computed in order to reproduce the lightest ++ and 
CT + states. The three states 2 ++ , ++ ,jmd 0~ + have been chosen because they are possible 
experimental glueball candidates [lsl. 



small errors in lattice calculations 



5] 



19] and because they are computed with relatively 



B. Confinement potential 



A good approximation of the confining interaction between a quark and an antiquark 
in a meson is given by the linear potential ar, where r is the distance between the two 
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particles and where a is the string tension. In a baryon, lattice calculations and some 
theoretical considerations indicate that each quark generates a flux tube and that these 
flux tubes meet in a junction point Rq which minimizes the potential energy (the so-called 
Y-junction). Following this hypothesis, the confinement in a baryon could be simulated by 
the three-body interaction 

3 

V qqq = a^2\vi- R \, (1) 

i=l 

For such a potential, the point Rq minimizes also the length of the three flux tubes and is 
identified with the Toricelli point 211 ] . 

The energy density A c of a flux tube (string tension) can depend on the color charge c 
which generates it. Lattice calculations 22j and effective models of QCD 23[ predict that 
the Casimir scaling hypothesis is well verified in QCD, that is to say that the energy density 
is proportional to the value of the quadratic Casimir operator F c 2 of the colour source 

A C = F C V (2) 



We have then A r , 



A,-, 



4er/3 = a and X g = 3 a. In this work we will assume that the 



confinement in a three-body color singlet is given by 



V, 
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E^ 2 



Rn 



(3) 



i=i 



This potential can be considered as the three-body generalization of the confinement used 
in Ref. [14]. No constant potential is added, contrary to usual Hamiltonians in mesons and 



baryons [7, 24]. Let us note that if the three color charges are not the same, Rq is no longer 
identified with the Toricelli point 251 ] . 

Interaction ([3]) is very difficult to use in a practical calculation. A good approximation 
can be obtained for three identical color charges by substituting Rq by the center of mass 
coordinate R cm and by renormalizing the potential by a factor / which depends on the 
three-body system J21J]. For three identical particles, the best value is f — 0.9515. We 
will use this approximation in the following, which seems more realistic than a confinement 
obtained by the sum of two-body forces |9j, lid ]. 

We include in our model the spin-orbit correction to the confinement potential. It is 
given by the Thomas precession of the particles and reads 261 ] 
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(4) 



This second order correction depends on the effective gluon mass fi. We review in Sec. Ill CI 
the main feature of this effective mass. Let us note that a Y-junction for the confinement 
leads to a three-body spin-orbit since the Torricelli point is function of the three particles. 
But with our approximation relying on the center-of-mass, the Thomas precession term is a 
sum of three one-body spin-orbit interactions. 

In Refs. 101 ]. the confinement potential saturates at large distances in order to 



simulate the breaking of the color flux tube between gluons due to color screening effects. 
An interaction of type seems a priori inappropriate since the potential energy can grow 
without limit. But the phenomenon of flux tube breaking must only contribute to the 
masses of the highest glueball states. Moreover, it has been shown that the introduction of 
a saturation could not be the best procedure to simulate the breaking of a string joining two 
colored objects [2j 



C. Dynamical constituent gluon mass 
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Within the auxiliary field formalism (also called einbein field formalism) 
be considered as an approximate way to handle semirelativistic Hamiltonians 
effective QCD Hamiltonian has a kinetic part depending on the particle current masses rrii 
and the interaction is dominated by the confinement. A state-dependent constituent 
mass /J,i = (>ypf + mf) can be defined for each particle, and all relativistic corrections (spin, 
momentum, . . . ) to the static potentials are then expanded in powers of l//ij. This approach 
has been used in Ref. [ijj] to build the two-gluon Hamiltonian. So, the same formalism will 
be applied also in this paper. 

Taking into account the considerations of Sec. Ill B\ the simplest generalization to a three- 
gluon system of the dominant part of the model III two-gluon Hamiltonian of Ref. [l^] is 

3 „ 3 



Ha 



i=i 



(5) 



i=l 



with the condition Y^h=i Pi = 0> since we work in the center of mass of the glueball. The 
gluons have vanishingcurrent masses and their color is such that (Ff) = 3. Contrary to 



„aflM. b , u, m 



been shown that it is an important ingredient to obtain correct two-gluon glueball spectra. 



5 



Using the technique of Ref. 29], it is possible to obtain an analytical approximate formula 
giving the glueball mass M and the constituent gluon mass /i (the three constituent gluon 
masses are the same since the wave function is completely symmetrized, see Sec. IIIII) 



M ^6/i with y/fa(N + 3). (6) 

N — 0,1, ... is the excitation number. With a value of the meson string tension a = 4 a/3 
around 0.2 GeV 2 , the smallest gluon constituent mass is around 650 MeV. It is then relevant 
to use an expansion in powers of l/^o- Such a value of the gluon mass is in agreement with 
the values used in Refs. . \lo\. but here the constituent mass is dynamically generated. 
Instead of using the auxiliary field formalism, it is possible to consider relativistic correc- 



tions which are expanded in powers of l/E^Pi) where E^pi) = \/ pf + mf (see for instance 
Ref. jsij]). But, this leads to very complicated non local potentials which are difficult to 
handle. 



D. Short-range potential 

The Hamiltonian Hq ([5]) gives the main features of the three-gluon glueball spectra, 
but the introduction of a short-range potential is necessary to achieve a detailed study. In 
Ref. 14], a one-gluon exchange (OGE) interaction between two gluons, coming from Ref. 



has been considered. It is not possible to use it directly for a three-gluon glueball because 
the color structure of the interaction is different. So, we use here the last version of a 



OGE interaction between two gluons developed specifically for three-gluon glueballs [9l. Il0|. 
Its explicit form, which is very similar to the form of the OGE interaction for two-gluon 
glueballs, is given below. 

This interaction contains a tensor part and a spin-orbit part. In our previous study of the 
low-lying spectrum \v\ . we neglected these terms since we worked on L = states. It has 
been shown that the tensor interaction between two gluons is small in two-gluon glueballs 



14j . Hence, we only add a two-body spin-orbit interaction between each pair of gluons. 

The OGE two-gluon potential has a priori a very serious flaw: depending on the spin state, 
the short-range singular part of the potential may be attractive and lead to a Hamiltonian 
unbounded from below [loj]. This problem is solved, as in Ref. jl4|, by giving a finite size 



to the gluon (see Sec. Ill Ep 
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The OGE potential depends on the gluon constituent mass. To determine it, we follow the 
procedure proposed in Ref. [14j. For a given set of quantum numbers {a}, the eigenstate \<p a ) 
of the Hamiltonian Hq is computed. With this state, a constituent gluon mass is computed 
A*q = {4 > o\y^\4 > otj- This value of \i a is then used in the complete Hamiltonian (see Sec. Ill F|) 
to compute its eigenstate with quantum numbers {a}. It is worth noting that, with this 
procedure, two states which differ only by the radial quantum number are not orthogonal 
since they are eigenstates of two different Hamiltonians which differ by the value of \i. It is 



shown in Ref. 
weak. 



16| that this problem is not serious, the overlap of these states being generally 



E. Gluon size 

In potential models, the gluon is considered as an effective degree of freedom with a 
constituent mass. Within this framework, it is natural to assume that a gluon is not a pure 
pointlike particle but an object dressed by a gluon and quark-antiquark pair cloud. Such an 



hypothesis for quarks .eads to very good resu.ts in meson |3l| and baryon M sectors. As 
in Ref. [14| , we assume here a Yukawa color charge density for the gluon 



1 e- u ^ 

4-7T7 2 u 



where 7 is the gluon size parameter. The interactions between gluons are then modified by 
this density, a bare potential being transformed into a dressed one. 

The main purpose of the gluon dressing is to remove all singularities in the short-range 
part of the interaction [32]. But, for consistency, the same regularization is applied to the 
confinement potential, although no singularity is present in this case. We think that the 
definition of a gluon size, which has a clear physical meaning, is preferable to the use of a 
smearing function only for potentials with singularity j^. flp|. 

A one-body potential, like \ri — R cm \, is dressed by a simple convolution over the density 
of the interacting gluon and the potential 

V{r)* = J dr'V{r')p{r -r'). (8) 

A dressed two-body potential, depending on \r{ — Vj\, is obtained by a double convolution. 
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This procedure is equivalent to the following calculation 



V(rY* = [ dr' V(r') T(r - r') with T(u) = 

J 8^7^ 



-u/y 



(9) 



Note that for the spatial dependencies of the spin-orbit interactions, which depend on the 
derivative of the first order potentials, the convolution must be performed before taking the 
derivative. 

For non vanishing value of 7, the value of the confinement potential at origin increases. 
This shifts the whole spectra to higher masses. Moreover, the strength of the Coulomb 
interaction is reduced. This also implies an increase of the glueball masses. 

F. Total Hamiltonian 

To obtain the total Hamiltonian for three-gluon glueballs which is the simplest general- 

[ I 

ization of the Hamiltonian for two-gluon glueballs from Ref. [14], we take the Hamiltonian 
Hq given by the relation (jSJ) and its spin-orbit correction (jlj) ; we add the OGE interactions 
coming from Ref. [h]] (without tensor part); and we dress all the potentials with the gluon 
color density (JTj). This gives the following Hamiltonian (Ff — 3, Fi ■ Fj = —3/2) 
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(10b) 

(10c) 
(lOd) 

(lOe) 
(lOf) 
(10g) 



where Y^=i P« = an d $y = Si + Sj. (3 — +1 (—1) for a gluon pair in color octet antisym- 
metrical (symmetrical) state. The constituent state-dependent gluon mass \i is computed in 
advance with a solution of the Hamiltonian Hq. 
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III. WAVE FUNCTIONS 



A gluon is a I{J P ) = 0(1~) color octet state. Two different three-gluon color singlet 
states exist 9J], which are completely symmetrical or completely antisymmetrical. The total 
isospin state of a glueball is an isosinglet and is completely symmetrical. Different total spin 
states are allowed with different symmetry properties. They are presented in Table [Til As 
gluons are bosons, the total wave function must be completely symmetrical. Its parity is 
the opposite of the spatial parity, and its C-parity is positive for color antisymmetrical state 
and negative for color symmetrical state. Let us note that a two-gluon glueball has always 
a positive C-parity. 



In our previous work [13], we mainly considered glueballs with the lowest masses. These 
states are characterized by a vanishing total orbital angular momentum L = and by 
a spatial wave function completely symmetrical with a positive parity. This immedia tely 
implies that the lowest glueballs are states with J PC equal to CT + , 1 , and 3 0, flQ| . 
No firm conclusion could be drawn for the 0~ + state because it appears both in two-gluon 
and three-gluon spectra. Our results for the 1 and 3 glueballs were in good agreement 
with lattice masses. We also presented the mass of the 2 even though it was found with 
a mass higher than the one predicted by the lattice calculations. We could understand this 
bad result as a first hint that possibly our approach could not be the best way to handle 
pure gauge spectra. 

We now extend our study to states with negative C-parity. Indeed, at least 3 gluons 
are contained into those glueballs. Lattice QCD recently computed masses of J + ~ glue- 
balls: 0+-(4780),l+-(2980),2 + -(4230), 3 + "(3600) in Ref. [5|, and l+-(2670), 3+-(3270), 
5 + ~(4110) in Ref. jf]]. The positive parity requests an odd angular momentum. 

In order to explain the high energy behavior of the proton-antiproton scattering, the 
existence of a trajectory carrying vacuum quantum numbers was postulated: The pomeron 

34j . The matching of the two-gluon spectrum and this trajectory is an important success 
of the theory 6|, |35|]. The odderon is the negative C-parity counterpart of the pomeron 

3. One often argnes that the odd spin gluebaUa he on its Rectory. The fftst two states 
(1 and 3 ) on its trajectory were already computed in Ref. [171 ] ■ In this paper, we also 
computed the 5 mass to check this hypothesis. Although lattice do not report any result, 
we will compare our mass to a Coulomb gauge approach for the odderon [371] ■ The oddballs 
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J masses for J = 1, 2, 3, 5, 7 as well as the 0~ + are computed within the Coulomb gauge 
approach in this reference. 

In order to reach a good accuracy, the trial spatial wave functions are expanded in large 
Gaussian function basis [38] . Recently, we found the Fourier transform of a Gaussian function 



and applied it to give the matrix elements for the semi-relativistic kinetic energy 39]. The 



interested readers may also find the matrix elements for spin dependent operators in this 



base in Ref. 40], where applications for three-body systems were presented. With more 
than 10 Gaussian functions for each color /isospin/sp in channels, we have checked that the 
numerical errors on masses presented are around or less than 1 MeV. 

Using the a value from our previous models A and B [y], [itJ , eigenvalues of the Hamil- 
tonian Hq ([5]) are presented in Table UTTl for various J PC quantum numbers. For each state, 
the corresponding constituent gluon mass /i is indicated. It is used to define the complete 
Hamiltonian H (TTTTj) . Let us note that, even if the Hamiltonian H does not contain spin- 
dependent potential, the spin of the wave function can strongly influence the mass through 
the symmetry of its spin-part. With our numerical procedure, the J +± glueballs are char- 
acterized by L > 1 spatial wave functions with L odd. In Table IIIIl one can see that the 
four lowest J + ~ glueballs are degenerate (L — 1, S — 1) as long as no spin-dependent term 
is included in the potentials. 

IV. RESULTS 

We present here the three-gluon glueball masses obtained with the complete Hamiltonian 



H (JTUJ) together with the two-gluon glueball masses computed in Ref. [14J (see Table [TV 



These masses are compared with results obtained by the lattice calculations of Ref. [ij 
This last work is an update of a previous study [4|. So, states not computed in Ref. [5|] 
but presented in Ref. (jl are also considered here. The 5 ±_ masses are also computed and 

n 

compared to the result of Ref. [6] . 

This work is an extension of our previous study [17J] in which we concluded that our 
1 and 3 states were in good agreement with lattice results. However, our 2 state 
is clearly higher than its lattice counterpart. The lattice results predict a 2 state at 
4010 MeV near the 1 and 3 states. With our Hamiltonian, a mass more than 1 GeV 
above is computed. It is unavoidable in our model, since a spin 2 function has a mixed 
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symmetry which implies a mixed symmetry for the space function and then a greater mass 
for the corresponding glueball, in agreement with the results of Refs. ^, [loj]. It has been 
checked that the spatial wave function of the 2 state is dominated by a configuration in 
which each internal variable is characterized by one unit of angular momentum. However 
in Ref. 371 ] . the 2 state lies between the 1 and 3 . Actually, they computed the 
spin-dependent corrections in perturbation with a wrong wave function. When the correct 



44)- 



wave function is used the 2 goes higher in agreement with our model 

Moreover in our model, as long as we do not include spin-orbit forces the + ~, 2 + ~, 
and 3 +_ glueballs are degenerate L = 1 states with the corresponding mixed symmetry 
for the spin S = 1 or 2 (see Table HT1) . We included spin-orbit interactions (coming from 
the one-gluon exchange and Thomas precession) to split the degeneracy. But firstly, the 
splitting between these states is not sufficient. And secondly, the hierarchy is not correct. 
Indeed, we checked that the two spin-orbit contributions canceled roughly each other. The 
one coming from the OGE decreases the masses by 100-200 MeV when the one coming 
from the confinement increases the masses by around the same amount. This fact is well 



known in baryons [45[ and hence it is not very surprising in a glueball potential model. Let 



us note that this cance 
two-gluon glueballs [ll|, 



lation between two spin-orbit contributions was already noted for 



141 ]. but the effect coming from OGE was always stronger. In lattice 



calculations however, the situation is strongly different: The splitting between l +_ and H 
is about 1.8 GeV! Such a splitting cannot be reproduced in our simple model. 

The values of parameters a and as are quite well determined with lattice calculations and 
Regge phenomenology. The situation is very different for the gluon size 7. In this work, it is 
assumed that 7 has the same value in two-gluon and three-gluon glueballs. We have checked 
that variations of the gluon size for three-gluon glueballs does not change the hierarchy of 
states, but the masses are globally shifted. As expected, the masses increase with 7. 

It is worth noting that the 1 H and 3 +_ three-gluon glueballs have masses similar to 
J = 2 and J = 3 two-gluon glueballs. Such low masses are difficult to explain. It is possible 
that these states are actually two-glueball bound states 4J]. Let us also note that l +_ and 
3 + ~ glueballs have low masses in the closed flux tube model 4f|. 

We have no firm explanation for such discrepancies. These problems could arise because 
the gluon has a constituent mass within our formalism. So, it possesses a spin as any massive 
particle, that is to say three states of polarization, and we use a basis where each state is 
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labeled by a couple (L,S). In the two-gluon sector, the (L,S)-basis leads to more states 
that the ones predicted by the lattice and also to J = 1 states forbidden by Yang's theorem. 
A solution to cancel the extra states is to deal with the so-called helicity formalism \^\ . 
Recently, we reviewed and applied this formalism to two-gluon glueballs [48]. We showed 
that a simple Cornell potential together with the helicity formalism reproduces the correct 
hierarchy given by the lattice. An instanton induced forced was also needed to raise the 
degeneracy between the scalar and pseudoscalar glueballs. In lattice calculations, the gluon 
is a massless particle with a definite helicity and then only two states of polarization. We 
suspect that the implementation of the helicity formalism for three-body systems would 
solve the hierarchy problem in the J + ~ sector but also for the 2 . This difficult task is 
out of the scope of this paper and we leave it for future work. Nevertheless, a attempt to 



generalization of the helicity formalism for three-body system is presented in Ref. 



49 



50]. 



It has been suggested that a three-body interaction can inverse the state ordering in the 
gluelump sector [511] • It is not clear that such a force acts also in glueballs and could solve the 
hierarchy problem. Actually, the dominant three-body force in glueball is the confinement, 
supplemented by its spin-orbit correction. In this work, as explained above, the confinement 
and its spin-orbit correction are simulated with a sum of one-body interactions. 

The odderon trajectory carries the quantum numbers J with odd J. This trajectory 



was investigated in Ref. 



371 ] within a Coulomb gauge Hamiltonian formalism. In this work, 



states are built with well defined (L, S) quantum numbers. The spectrum is in good agree- 
ment with our results (see Table IV|) . Let us note that the Coulomb gauge Hamiltonian is 
more complicated than ours. For instance, the annihilation diagram for two gluons is taken 
into account whereas it vanishes in a potential model 52J. It is not surprising that both 
approaches gives rather the same results since they used the same formalism (a (L, 5*)-basis). 
We found that the two first oddballs 1 and 3 lie in lattice error bars. This fact can be 
surprising at first glance. Indeed, we invoked the helicity formalism as a possibility to solve 



the mass problem of the J + . Hence the same problem should arise for the J 



helicity formalism, a given 

jPC 

is a particular combination of (L,S) couples [43, l48| . We 



n the 



suspect than the oddballs are largely dominated by the component (L — J — 3, S — 3). It 
would be interesting to have more information about those states but unfortunately lattice 
have difficulties to identify higher excited states. 
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V. CONCLUSION 



The masses of pure three-gluon glueballs have been studied with the generalization of a 
semirelativistic potential model for pure two-gluon glueballs [jjj]. The short-range part of 
the potential is the sum of two-body OGE interactions. For the confinement, a potential 
simulating a genuine Y-junction is used and it is assumed that the Casimir scaling hypothesis 
is well verified. The gluon is massless but the OGE interaction is expressed in terms of a 
state-de pen dent constituent mass. The Hamiltonian depends only on 3 parameters fixed 
in Ref. [141 ]: the strong coupling constant as, the string tension a, and a gluon size 7. 



All masses have 



39 



)een accurately computed with an expansion of trial states in Gaussian 



functions [38|, 

/,From our previous paper 



4Q|. 

14, EtI . we know that J ±+ two-gluon glueball spectra is in good 
agreement with lattice calculations % Q , but extra states not seen in lattice calculations are 
predicted. The J _± three-gluon glueball are also in quite good agreement with these lattice 
results, except for the 2 state which is computed with a very high mass in our model. 
In this work, the 5 and some three-gluon glueballs are finally computed. The J 
candidates with J odd for the odderon trajectory are in agreement with some other works. 
But, if the +_ and 2 + ~ states are predicted in quite good agreement with lattices results, 
it is not the case for the l +_ and 3 H states computed with too high masses in our model. 
One could interpret such a discrepancy as due to the fact that the spin-orbit forces are too 
feeble in our model to raise the degeneracy between J + ~ states. But we do not believe that a 
physical process, ignored here, is able to produce a strong enough spin-orbit force giving rise 
to more than 1 GeV energy gap. We think that the strong discrepancies between our results 
and lattice computations are due to the fact that gluons do have constituent non vanishing 
masses in our approach. They are then characterized by a spin, and not by a helicity as it 
could be expected for particles with a vanishing current mass. We have shown that working 
with the helicity formalism could cure the problem of extra states in two-gluon glueballs 
481 ]. We think that this formalism applied to three-gluon glueballs could also improve the 
predictions of a potential model. Such a work is in progress. 
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TABLE I: Parameters for models A and B (a = 3 a/4). For both models, the gluon current mass 
is zero and / = 0.9515. 



Model A Model B 



a 0.16 GeV 2 0.21 GeV 2 

a s 0.40 0.50 

7 0.504 GeV" 1 0.495 GeV" 1 



TABLE II: Characteristics of three-gluon spin functions with total spin S, intermediate couplings 
S'int, and symmetry properties which can be obtained by coupling (A: Antisymmetrical, S: Sym- 
metrical, MS: Mixed symmetry). The multiplicity of each symmetry type is indicated. 



S Sint Symmetry 

1 1 A 

1 0, 1, 2 1 S, 2 MS 

2 1, 2 2 MS 

3 2 IS 
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TABLE III: Masses Mq of the Hamiltonian Ho (0) as a function of the J PC quantum numbers. 
The (L, S) quantum numbers are indicated with the corresponding constituent gluon masses n$. 
Values in MeV are computed with the value of a from models A/B. The lowest masses are printed 
in italic. 



J PC 


(L,S) 


M 


Mo 


J PC 


(L,S) 


M 


Mo 


n — 
u 


IU,U ) 


DO 1 4/ DooD 


QOQ / 1 nfi/l 

y/y/ 1UD4 


u 


(0,0J 




CQC /CI Q 

OoD/Dlo 


1~ 


(0,1) 


3211/3679 


535/613 


1-+ 


(0,1) 


4156/4761 


693/794 


2~ 


(0,2) 


4156/4761 


693/794 


2-+ 


(0,2) 


4156/4761 


693/794 


3~ 


(0,3) 


3211/3679 


535/613 


3-+ 


(0,3) 


5574/6385 


929/1064 


5~ 


(2,3) 


4182/4791 


697/795 










o+- 


(1,1) 


3752/4298 


625/717 










1+- 


(1,1) 


3752/4298 


625/717 










2+- 


(1,1) 


3752/4298 


625/717 










3+- 


(1,2) 


3752/4298 


625/717 










5+- 


(3,2) 


4596/5265 


766/878 
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TABLE IV: Glueball masses in MeV and (glueball mass ratios normalized to lightest 2 ++ ). The 
two-gluon masses are taken from Ref. The error bars for lattice mass ratios are computed 

without the normalization error on the masses. The lightest ++ , 2 ++ , and ^ states are taken 
as inputs to fix the parameters. The first column indicates the valence gluon content as predicted 
by our model. 



J PC Lattice [Ref.] Model A Model B 



gg o++ 


1710±50±80 


(0.72 ±0.03) 


M 


1604 


(0.78) 


1855 


(0.78) 




2670±180±130 


(1.12 ±0.09) 


[4] 


2592 


(1.26) 


2992 


(1.26) 


2++ 


2390±30±120 


(1.00 ±0.03) 


[5] 

L J 


2051 


(1.00) 


2384 


(1.00) 


o-+ 


2560±35±120 


(1.07 ±0.03) 


[51 

1 J 


2172 


(1.06) 


2492 


(1.05) 




3640±60±180 


(1.52 ±0.04) 


1 J 


3228 


(1.57) 


3714 


(1.56) 


2-+ 


3040±40±150 


(1.27 ±0.03) 


[5] 

L J 


2573 


(1.25) 


2984 


(1.25) 




3890±40±190 


(1.63 ±0.04) 


[41 

1 J 


3345 


(1.63) 


3862 


(1.62) 


3++ 


3670±50±180 


(1.54 ±0.04) 


[5] 

1 J 


3132 


(1.53) 


3611 


(1.51) 


ggg i 


3830±40±190 


(1.60 ±0.04) 


[5] 

L J 


3433 


(1.67) 


3999 


(1.68) 


2— 


4010±45±200 


(1.68 ± 0.04) 


[5] 

1 J 


4422 


(2.16) 


5133 


(2.15) 


3~ 


4200±45±200 


(1.76 ±0.04) 


[5] 


3569 


(1.74) 


4167 


(1.75) 


o-+ 








3688 


(1.80) 


4325 


(1.81) 


o+- 


4780±60±230 


(2.00 ± 0.05) 


[5] 


4043 


(1.97) 


4656 


(1.95) 


1+- 


2980±30±140 


(1.25 ±0.03) 


[5] 


3992 


(1.95) 


4626 


(1.94) 


2+- 


4230±50±200 


(1.78 ±0.04) 


[5] 


3907 


(1.90) 


4542 


(1.91) 


3+- 


3600±40±170 


(1.51 ±0.04) 


[5] 


4033 


(1.97) 


4568 


(1.92) 


5+- 


4110±170±190 


[6] 




4571 


(2.23) 


5317 


(2.23) 


5~ 








4521 


(2.20) 


5263 


(2.21) 
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TABLE V: Odderon quantum numbers and masses in MeV. L and S assignments agree in our 
model B and in Ref. [371 ] . 





1— 


3— 


5~ 


7— 


s 


1 


3 


3 


3 


L 








2 


4 


Model B 


3999 


4167 


5263 




Coulomb Gauge [37] 


3950 


4150 


5050 


5900 


Lattice [5] 


3830 


4200 






Lattice [6] 


3100 


4150 






Wilson loops [53] 


3490 


4030 







20 



> 4000 



Lattice (2006) 
Lattice (2006) - refercne 
Lattice (1999) 
Model B 



•I 9 



2 2 3 



o 



9 99 



\ 2 .1 1 2 3 



FIG. 1: Glueball masses given in MeV. Dotted diamonds: Results from model B (two-gluon masses 



n 



are taken from Ref. [14[] and three-gluon J masses are taken from Ref. [13]); Black and white 
circles: Lattice results from Ref. [5(; White squares: Uttc^^U. Black circles 
indicate the reference states taken as inputs to fix the parameters. The error bars for lattice results 
are computed by summing the two uncertainties (see Table HVj) . 
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